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ABSTRACT 

In this manuscript, we study properties of BLR of well-known reverberation-mapped AGN, in order 
to find reliable evidence for intermediate BLR. We firstly check properties of mapped AGN collected 
from literature in plane of ^\pl<^\ a vs ^blr/^blr- Commonly, virial BH masses based on observed 
broad Ha and H/3 should be coincident. However, among the mapped objects, PG0052 and NGC4253 
are two apparent outliers in the plane of c^/a^ vs 

i?|g R /i?BLR> which indicate BLRs of PG0052 
and NGC4253 have some special characters. Then based on the 55 public spectra of PG0052, BLR of 
PG0052 is been carefully studied in detail. We find that line width ratio of total observed broad Ha to 
total observed broad H/3 is ~0.7, which is much smaller than theoretical/observational value of ~0.9. 
Furthermore, flux ratio of total broad Ha to total broad H/3 is about 6.8 (Balmer Decrement), which is 
not one reasonable value for BLUE quasar PG 0052+251. Moreover, properties of line cores based on 
PCA technique confirm there is one inner broad component and one seriously obscured intermediate 
broad component in BLR of PG0052. If the seriously obscured intermediate BLR was accepted, prop- 
erties of PG0052 in the plane of o\plo\ a vs -Rblr/^blr could be well reproduced, which indicates 
that the intermediate BLR actually is well appropriate to mapped quasar PG 0052+251. Finally, the 
large distance between inner component of BLR and intermediate component of BLR based on CCF 
results rejects the possibility that the intermediate component is probably extended part of inner 
component of BLR. 

Subject headings: Galaxies:Active - Galaxies: nuclei - Galaxies:quasars:Emission lines - Galaxies: 
Individual: PG 0052+251. 



1. INTRODUCTION 

Through strong broad emission lines coming from 
Broad Line Regions (BLR) of Active Galactic Nu- 
clei (AGN), properties of BLR of AGN which can 
not be directly resolved by current observational 
technique and instruments have bee n understood 
deepe r and deepe r (see review s of ISulentic et al 
2000: IGaskelll 1201(1 iKrause et al l I20TD: IDennev et al 



[20061 [2001 IDennev et all [20Toh based on the rever- 
beration mappi ng technique (|Blandford fe Mckeel 119821 : 
IPetersonl I1993D . There are so far more than 40 
nearby broad line objects, of wh ich .Bbt.r and virial 



BH masses have been determined (Peterson et al 
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2004 



Dennev et alJl201tt iKaspi et alJl2005t iBentz etaT.112009 
2010t iBarth et aLT 20111 ). Then one empirical relation 



20101: iNetzer fc Marzianil 120101: iSluse et all 120111 : 
Pancoast et al.ll2011l etc.. and references therein). Based 
on properties of BLR of AGN (especially virialized 
emission line clouds in BLR.lGaskel]||1988trWandel et all 
19991 : iPeterson fc Wandell fl999t lOsterbrock fc Mathews 
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19861 : IDennev et al.ll2009T lGaskcll 200j), the most conve- 



nient method to es timate virial BH masses of broad lin e 
AGN is proposed dVestergaard 12001 lOnken et al.ll200 
Peterson et al.1 120041: IPetersonl 120101: IPeterson fc Bent: 
20061: ICollin et al.ll2006TTKellv fc Bechtol d 2007) 



M BH «r X i?BLR 



(1) 



where V represents probale rotating velocity of emis- 
sion line clouds in BLR (line width of broad emission 
line), i?BLR means distance between BLR and central 
black hole of AGN (size of BLR) which can be mea- 
sured by time lag between broad line emission and 
contin uum emission from long-period ob s erved spec- 
tra (jKaspi et al.ll2000l: IPeterson et alj 120041 : IBentz et all 
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has been found for -Rblr, -Rblr tx L 
L5100 means AGN continuum lumino sity at 5100Aand 
Lun e . represen ts broad line lumino sity dKaspi et al.ll2005l: 
Dennev et all 120101: IBentz et all l2009t iWang fc Zhand 
20031 : iGreene et al.ll2010f ). and virial BH masses can be 
simply and conveniently estimated from single epoch 
spectra of broad line AGN by line width and contin- 
uum l uminosity (or broad line luminosity) dVestergaarc ' 
2001 iMcLure fc Dunloijl [200l ISulentic et all 1 200? 
Netzer fc Marzianil 12010: ISulentic et alj 120061 12007 
Shen et al.ll2008t iRafiee fc Halll20lH INetzer k Marzianil 
2010t iKellv fc Bechtoldl 120071 : IGreene fc Hoi 120051 Den- 
ney et al. 2009a) 

M BH cx V 2 x L~ 00 5 ex V 2 x L'ifJ (2) 

It is clear that from single epoch spectra of broad line 
AGN, virial BH masses based on different broad emis- 
sion lines should be coincident. However, some conflict - 
ing results have been reported. Kelly & Bechtold (2007) 
have found that virial BH masses based on broad lines in 
UV band and broad lines in optical band are coincident 
well (so me simple and more recent resul ts can also be 
found in IRafiee fc Halll |20TU ) . However. ISulentic et al.1 
(|2007f) have shown that UV bro ad line CIVA1545 A is a 
poor virial estimator. Recently, Assef et al. ( 2010) have 
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found the opposite that virial BH mass estimates Based 
on CIVA1545A are consistent with those based on the 
balmer lines. Different virial BH masses based on differ- 
ent broad emission lines probably provide some kindly 
deep information of structures of BLR of AGN. In this 
manuscript, to study properties of virial BH masses on 
different broad emission lines for reverberation-mapped 
AGN is our main objective. In order to find more reliable 
results, parameters based on the reverberation-mapped 
AGN rather than AGN with single epoch spectra should 
be firstly considered. For recent reported reverberation- 
mapped AGN, sizes of BLR are based on long-period 
varied broad line emission lines in optical band, espe- 
cially balmer lines. Thus, in the manuscript, we mainly 
consider the properties of balmer lines (Ha and H/3) of 
reverberation mapped AGN. 

Besides the reverberation mapping technique and viri- 
alization method for virial BH masses, there are some 
references which study properties of e mission lines re- 
gions for broad optical balmer lines. iKorista fc Goadl 
(2004) have shown that after consideration of luminos- 
ity dependent responsivities, theoretically estimated size 
of BLR based on Ha should be 20% larger than size 
of BLR based on H/3, and corresponding line width of 
H/3 should be some larger than that of Ha. The the- 
oretical results can be confirmed by observational re- 
sults shown in iBentz et al.l (|2010l ) (size of BLR on Ha 
is some large r than size of BLR o n H/3) for mapped ob- 
jects and in iGreene fc Hoi (|200l (1 ine width of broad 
H/3 is some larger than line width of broad Ha) for pure 
quasars. Furthermore, one cal led two-component model 
of BLR has be en proposed (iPopovid 120071 : iBon et all 
120091: iZhu et al.ll2009t IHu et alJl2008[ ): one inner region 
for broad wings and one intermediate region for cores, 
in order to explain complexed broad lines of AGN. Al- 
though estimated (or measured) parameters (line width 
and size of BLR) of Ha and H/3 are some different to some 
extent, estimated virial BH masses based on Ha and H/3 
sh ould be coincident, s uch as the shown confirmed results 
in IGreene fc Hoi (|2005l ) for pure quasars. 

It is interesting to check the virial BH masses based on 
Ha and H/3 for reported reverberation mapped AGN. In 
the manuscript, we will find that there are some mapped 
AGN, virial BH masses based on Ha are much different 
from BH masses on H/3, which will indicate some spe- 
cial information for structure of emission line regions for 
balmer lines. This manuscript is organized as follows. 
Section 2 shows the results based on parameters (mea- 
sured line width and size of BLR based on long-period 
variations of line emission and continuum emission) of re- 
ported reverberation mapped AGN. Section 3 gives the 
detailed results for blue quasar PG 0052+251 with much 
different virial BH masses on Ha and on H/3. Section 4 
shows the discussion and conclusion. In this manuscript, 
the cosmological parameters Hq — 70km • s _1 Mpc _1 , 
SIa = 0.7 and Vt m = 0.3 have been adopted here. 

2. VIRIAL BH MASSES ON BALMER LINES FOR 
MAPPED AGN 

As shown in Introduction, we try to check 
virial BH masses on Ha and H/3 for mapped 
AGN. The reported mapped AGN are mainly 
included in two research groups, AGNWATCH 



group (public data of 10 AGN can be foun d in 
http: / / www, astronomy, ohio-state . edu / ~ agn watch /) 
and research group in Wise Observatory in Tel Aviv 
University (public data of 17 PG quasars c an be found 
in |http:/ /wise-obs. tau.ac.il/~s hai/PG/, Kas pi et al.1 
2000). Corresponding references about observational 
techniques and instruments for the two groups can be 
found in the two websites above. 

Before proceeding further, we rou ghly check results 
reported in the literature (especially iKaspi et al.l I2000L 
120051: iPeterson et~aT1 120041: IBentz et al.l 120101 ) for all the 
43 reverberation mapped AGN, of which observation 
data are public or unpublic. Here we mainly com- 
pare properties of broad Ha and broad H/3 (the two 
most strongest broad emission lines in optical band of 
AGN), including line widths (second moment rather than 
FWHM, which is the best approximation o f V in Equa- 
tion 1 as supposed by IPeterson et al.l 120041 ) and sizes of 
BLR (measured time lag between line emission and con- 
tinuum emission, not the calculated value from empirical 
relation of R oc L~ 5 ) based on Ha (^blr) and on H /3 
(^blr)i wn i cn are collected and listed in Table 1. Here, 
only objects with reliable parameters (line width and size 
of BLR), P > 1.5 x P orr (where P means measured pa- 
rameter, P crr represents corresponding uncertainty for 
the parameter) are collected from the literature. Based 
on the simple criterion, some reverberation-mapped ob- 
jects should be rejected. For example NGC3227, its size 
of BLR based on H/3 is R^ R ~ 8.2±| flight - days, 
uncertainty of 8.4 is larger than measured value 8.2, 
thus NGC3227 is rejected. We also should note that 
the famous reverberation-mapped object NGC5548 is re- 
jected due to t h e follo wing reasons. On the one hand, in 
IPeterson et al.l (|2004|) , line width of broad Ha is not re- 
liable due to larger uncertainty. On the other hand, we 
can find that size of BLR are similar, but line width of 
broad H/3 ( both second moment a nd FWHM) are much 
different in IPete r son et all (|2004f ) from those shown in 
IBentz et al.l (12010ft. ~ 2 000km/s and FWHMCH/31 ~ 
5800km/s in lPeterson e t al.l (12004ft. g Hft ~ 4200km/s and 
FWHM(H^) ~ 12000km/s in IBentz et all (|2010l) . Thus, 
NGC 5548 is not included in our parent sample listed in 
Table 1. Furthermore, there is another object, PG0844, 
we should note. For PG0844, there are reliable param- 
eters of line width and s ize of BLR for broad Ha and 
broad H/3 in lKaspi et al.l ()2000f ). however, there i s no re - 
liable size of BLR based on H/3 in lPeterson et al.l (|2004l) . 
Thus PG 0844 is also rejected. Eventually, there are 16 
objects listed in Table 1, 9 objects have public spectra 
(7 PG quasars observed by Wise Observatory and 2 ob- 
jects included in AGNWATCH project), and the other 7 
objects have no public spectra. Furthermore, we should 
note listed values of parameters in Table 1 are values from 
more recent literature. For example, values for PG0052 
can be found in Kaspi et al. (2000) and in Peterson et al. 
(2004). Then the listed values for PG0052 are collected 
from Peterson et al. (2004) . 

As shown in introduction, it is clear that there is one 
strong correlation for mapped AGN based on virial BH 
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masses 

2 v nH« _ 2 v r>Hf3 
a Ha x -"-BLR — °H/3 x ^BLR. 

R m (3) 
; "H« n2 _ -"-BLR v ' 

In this section, we will check the correlation for the 
selected mapped AGN. Figure [1] shows the correlation 
between line width ratio of broad H/3 to broad Ha 

(oh^/ojjq,) and size ratio of -Rblr to f° r tnc 16 

reverberation-mapped objects listed in Table 1. For all 
the 16 objects, the Spearman Rank correlation coeffi- 
cient for the correlation is about 0.4 with P nu \\ ~ 14% 
(two-sided significance of deviation from zero) , which in- 
dicates there is one rough positive correlation. In the 
figure, we also show the corresponding 99% confidence 
bands for the linear correlation u\^ja\ a = -Rblr/^blr- 
Based on the results in FigureHJ there are probable 2 out- 
liers, PG0052 and NGC4253, which are deviating from 
the linear correlation supported by virialization method. 
After the two outliers are rejected, the re-calculated 
Spearman Rank correlation coefficient is about 0.9 with 
P null - 4 x 10- 5 . 

Based on the results shown in Figure [TJ the much dif- 
ferent virial BH masses based on Ha and H/3 should in- 
dicate that there are some unique characters for broad 
balmer emission line regions. To find and study the 
character is the main objective of our following sec- 
tion. For the two outliers, only PG0052 have pub- 
lic observed spectra which can be downloaded from 
|http:/ / wise-obs.tau.ac.il/~shai/ PG/. This, in the fol- 
lowing section, we will mainly study the properties of 
reverberation mapped AGN PG0052, and try to find 
some special characters of broad balmer line regions of 
PG0052. 

3. MAIN RESULTS FOR PG0052 

In this paper, all the observation al data and s pectra 
of PG 0052+251 arc collected from Kaspi et al. (2000) 
( [http://wise-obs.tau.ac.il/~shai/ PG/). The detailed de- 
scription about th e data and spectra can be found in 
IKaspi et all (pOOO) . Here, we do not describe instru- 
ments and observational techniques any more. There 
are 56 spectra from 4000A to 8000A with dispersion of 
~3.8A/pixel and spectral resolution about ~ 10A ob- 
served from 16th Oct. 1991 to 27th Sep. 1998 for PG 
0052+251. However, there is one spectrum with many 
bad pixels around Ha, thus only 55 spectra are consid- 
ered. The collected spectra from the website above have 
been binned into lA per pixel and been padded from 
3000A to 9000A as shown in the website. We mainly 
consider the data and spectra as follows. 

3.1. Properties of Mean Spectrum 

Mean spectrum of PG 0052+251 can be created by 
PCA (Principal Components Analysis or Karhunen- 
Loeve Transform method, ) technique applied for ob- 
served noisy spectra. PCA technique is a mathemat- 
ical procedure that uses an orthogonal transformation 
to convert a set of observations of possibly correlated 
variables into a set of values of uncorrelated variables 
called principal components. Certainly, mean subtrac- 
tion (or mean centering) is necessary for performing 



PCA to ensure that the first principal component de- 
scribes the direction of maximum variance. However, if 
mean subtraction is not performed, the first eigencompo- 
nent through PCA technique commonly represents the 
mean spectrum of noisy spectra. Here, the convenient 
and public IDL PCA program 'pca_solve.pro' written 
by D. Schlegel in Princeton University is used, which 
is included in SPSS software package of IDLSPEC2D 
( |http: / /spectro.princeton.edu/[ ) . 

Figure [2] shows the mean spectrum of PG 0052+251 
with relative flux density flux(5100A) — 1. Further- 
more, the mean s pectrum of PG 0052+251 created by 
IKaspi et al.l (|2000| ) is also shown in our Figure [2j which 
is the same as our mean spectrum created by PCA 
method. Then balmer emission lines arc fitted by simple 
gaussian functions through Levenberg-Marquardt least- 
squares minimization method. The best fitted results are 
also shown in Figure [2j Emission lines around H/3 can be 
best fitted by one broad gaussian function for broad H/3 
and three narrow gaussian functions for narrow H/3 and 
[O m]A4959,5007A doublet. Here when [O lll]doublet 
is fitted, we require that [O iii]A4959, 5007A have the 
same line width in unit of km/s, and flux ratio of 
[O IH1A4959A to \Q IH1A 5007A is theoretical value 0.33 
(|Dimitriievic et aLll2007l ). Emission lines around Ha can 
be fitted by one broad gaussian function. Here we do 
not consider narrow emission lines around Ha in the 
mean spectrum, due to much weak [O i]A6300,6363A, 
[N n]A6548,6583A and [S n]A6716,673lA doublets. Fur- 
thermore, we should note that features of atmospheric 
A band near 7620A in observed frame can be detected, 
which are marked in Figure [2] and in the following Fig- 
ure[3]by symbol ffi. Here, we simply discuss effects of the 
features on measured line parameters. Line parameters 
around Ha are firstly measured without any considera- 
tion of effects of the features (functions are applied to 
total observed line profile), and then re- measured with 
consideration of effects of the features (functions are ap- 
plied to observed line profile with features of atmospheric 
A band rejected). Measured line parameters due to the 
two different procedures are similar. Thus, in the follow- 
ing part, effects of features of atmospheric A band near 
7620A are totally ignored. The measured line parame- 
ters of broad Ha and broad H/3 in mean spectrum are 
listed in Table 2. 

Based on the line parameters listed in Table 2, it 
is clear that line width (second moment) of broad H/3 
(<xh/3 ~ 38A ~ 2370km/s) is much different from line 
width of broad Ha (cthq ~ 37A ~ 1700 km/s), which 
is muc h differ ent from the results sho wn in iGreene fc Hoi 
([2005T) and in lKorista fc Goad! (|200l (a^/o-Ha ~ 1.06 — 
1.1). The results indicate that broad balmr line regions 
of PG0052 probably have some special characters which 
are some different from common quasars. In other words, 
single broad gaussian function applied to fit broad Ha 
and H/3 is not so good. Thus, we consider broad balmer 
emission lines by a different method as follows. 

3.2. Properties of Line Cores 

In order to more clearly show properties of broad 
balmer emission lines, PC A technique is applied a gain 
as what have been done in iBrotherton et all (|1994[ ) and 
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in iFrancis et al.l ((1992). Before to study properties of 
line cores, the commonly accepted step of mean subtrac- 
tion (or mean centering) is firstly performed. Then af- 
ter the PCA technique applied to the spectra with zero 
mean, the first principal component represents emissio n 
line cores (jBrotherton et all 119941 : IFrancis et al.l fl992) . 
Figure [3] shows emission line cores around H/3 and Ha. 
In the figure, line profile of H/3 with flux density scaled 
by 3 (intrinsic flux ratio of Ha to H/3 for blue quasar) 
is directly compared with line profile of Ha in the first 
principal component. Based on the results shown in Fig- 
ure [3j it is clear that the line cores of H/3 and Ha are 
much different, there is one apparent broad component 
in Ha but no broad component in H/3 in the first prin- 
cipal component. The results indicate besides the sim- 
ilar inner broad component as that of observed broad 
H/3, there is one other intermediate broad component in 
Ha but no intermediate broad component in H/3 (one 
seriously obscured intermediate component for balmcr 
lines). Moreover, we should note that there are narrow 
lines around H/3 in the first principal component shown 
in Figure [3J which is due to much extended component 
of [OIII] A4959, 5007A Similar narrow lines can also be 
found in rms spectrum of PG 0844 shown in Figure 2 in 
Kaspi et al. (2000). Those narrow lines can not affect 
our results about line cores of Ha and H/3. 

Thus, broad Ha is fitted again by two gaussian func- 
tions, one broad gaussian function with similar line 
width as that of broad H/3 and one intermediate broad 
gaussian function, through Levenberg-Marquardt least- 
squares minimization method. The two components are 
also shown in bottom-right panel in Figure [5] Param- 
eters of the two components of Ha are also listed in 
Table 2. The line width (second moment) of interme- 
diate broad component of Ha is about llOOkm/s which 
is much larger than the line width (second moment) of 
narrow [O iii]A5007 (<r([0 Hl])~ 6.6A=395km/s). Based 
on flux ratio of inner broad Ha to intermediate broad Ha, 
~ 2.6 — 3, it is clear that intermediate broad component 
of Ha is distinct and clearly decomposed in the mean 
spectrum of PG 0052+251. Certainly, we also try to fit 
broad H/3 by two broad gaussian functions. However, 
intermediate broad component of H/3 is not reliable due 
to smaller measured line width and line luminosity than 
corresponding measured uncertainties for the expected 
intermediate component of H/3, which is consistent with 
results shown in Figure [3l The results indicate the inter- 
mediate broad component of BLR probably exists and is 
seriously obscured (no intermediate broad component of 
H/3) for PG 0052+251. Certainly, the intermediate broad 
component of H/3 can not be detected perhaps due to low 
quality of spectra for PG 0052+251. Although there is 
no true value of flux ratio of intermediate broad compo- 
nent of Ha to intermediate broad component of H/3, it is 
not difficult to determine cut-off value of S/N (signal to 
noise) for future high quality of spectra. If we accept flux 
ratio of intermediate broad component of Ha to interme- 
diate broad component of H/3 is about 6, S/N should be 
larger than 24. 

Before the end of the subsection, we discuss proper- 
ties of %2 for fitted results for Ha and H/3, in order to 
find more reliable evidence for the existence of the second 
gaussian component for Ha. Here %2 is value of summed 



squared residuals for measured parameters divided by de- 
gree of freedom, which can be used as one good indicator 
to determine whether fitted results are acceptable. Here, 
we accept uncertainty for flux is about ~ 10%. The val- 
ues of x 2 are 0.18, 0.06 and 0.08 for best fitted results 
by one gaussian functions and by two gaussian functions 
for Ha and by gaussian functions for lines around H/3 
respectively. Then based on fitted results by one or two 
gaussian functions for Ha, F-test is performed: at 99% 
confidence level, calculated variance ratio of results by 
one gaussian function to results by two gaussian func- 
tions is 2.14 for Ha in mean spectrum, which is much 
larger than F value 1.34. The results indicate two gaus- 
sian functions for Ha is necessary. Due to the unreliable 
second gaussian component for H/3 through least-squares 
minimization method, F-test is not performed for results 
of H/3. 

3.3. Properties of Observed Spectra 

We cons i der observed spectra collected from 
iKaspi et al.1 ((2000) in the subsection. Line param- 
eters are measured through gaussian functions for all 
the 55 observed spectra with both Ha and H/3 for PG 
0052+251, as what we have done above for balmcr 
lines in the mean spectrum. Emission lines around H/3 
are measured once: only a single broad component is 
fitted for H/3. Emission line Ha is measured twice, 1): 
one broad gaussian function is fitted for Ha, 2): two 
broad gaussian functions are fitted for Ha, one broad 
component and the other intermediate broad component 
as shown in subsection above. When two gaussian 
functions are applied to fit Ha, broad component has 
similar line width as the one of broad H/3 (with a 
permitted scatter of O.ldex) and line width ratio of 
the corresponding two functions (broad to intermediate 

broad) is fixed to 2.16 (~ nookm/s ' based on the 
measured line widths of two components of Ha in the 
mean spectrum). 

After the measurements of line parameters, we can 
find if one gaussian function is applied to fit Ha, line 
width ratio of total broad H/3 to total broad Ha is about 
1.4 + 0.1, which is m uch larger than the m ean value ~1.1 
for quasars shown in lGreene fc Hoi (|2005[) and larger than 
the th eoretical value ~ 1.06 — 1.1 in IKorista fc Goadl 
(2004). Furthermore, if one gaussian function is applied 
to fit Ha, flux ratio of total broad Ha to total broad 
H/3 (Balmer decrement) is about 6.8, which is a larger 
and unreasonable value for blue quasar PG 0052+251 
(mean value from composite spe ctra of quasar is about 
3.56 in lVanden Berk et al.l l2001). However, if two gaus- 
sian functions are applied to fit Ha, flux ratio of inner 
broad Ha to inner broad H/3 (observed total broad H/3 
only includes the inner broad H/3, intermediate broad 
H/3 is seriously obscured) is about 3.61 which is consis- 
tent with the value from composite spectra of quasars 
(jVanden Berk et aT]|200lD . and line width ratio is about 
fHa/fH/3 ~ 0.94 (here Ha is the one without intermedi- 
at e broad Ha) , which is consist ent with the resul t show n 
in IGreene fc Hoi ((20051) and in IKorista fc Goadl (l200l . 
H/3 is slightly larger than Ha ((Thq/^h/3 ~ 0.91). Cer- 
tainly, we should note that there are some cases that 
intermediate broad Ha is too weak to be detected in ob- 
served spectra, as one example shown in Figure [4j which 
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is probably due to the observed spectra with lower res- 
olution. It is clear that total observed broad Ha sepa- 
rated into two components (one inner broad component 
and one intermediate broad component) should be more 
reasonable for blue quasar PG 0052+251. 

3.4. CCF Results 

Finally, we consider results from cross correlation func- 
tion (CCF function) applied to measure size of BLR of 
PG 0052 +251, which can be determined by time lag be- 
tween variations of continuum emission and variations of 
broad lines emission. Here the variations of broad emis- 
sion lines have four components, variations of observed 
broad H/3 (inner broad H/3), variations of observed total 
broad Ha, variations of inner broad Ha (coming from in- 
ner BLR) and variations of intermediate broad Ha (com- 
ing from intermediate BLR). Here the used flux densi- 
ties of total b road Ha, H/3 and c ontinuum emission are 
collected from iKaspi et al.l (|2000l ). And then, according 
to the measured flux ratio of broad Ha to intermediate 
broad Ha of each observed spectrum, the CORRECTED 
flux density of total broad Ha for each observed spectrum 
can be separated into two values for inner broad compo- 
nent and for intermediate broad component of Ha. Thus, 
effects from different instru ments in different configura- 
tions can be totally ignored (|van Groningen fc Wanders! 
1992). In other words, the following used flux densities 
of broad line components are reliable. 

Here, co mmon interpolated cross-correlation func- 
tion JI CCF) (IGaskell fc SparkdlroM IGaskell fc Peterson! 
Il987t iPetersonl |1993|) is applied to quantify time lag 
between continuum emission and broad lines emission. 
We do not co nsider z-transfer discrete correlation func- 
tion (ZDCF. IMexanderl fl997t lEdelson fc Krolikl Il98l 
IWhite fc Petersonlll994D any more, because results from 
ZDCF are exce l lent a g reeme n t with results from ICCF 
(IPeterson et al.l [l99l H992l [2001 IKaspi et al.l 120001 
iBentz et al.li201fj|r Figure [5] shows the final results based 



on the four broad components of Ha and H/3, includ- 
ing measured sizes of BLR and corresp onding uncertain- 
ties determined by b ootstrap method (jPress et al.lll992t 
IPeterson et al.lfl998h . Measured sizes of BLR based on 
variations of different broad components are listed in Ta- 
ble 2. The size of inner BLR is about 110 light-days based 
on inner broad Ha (or observed broad H/3) and the size 
of intermediate BLR based on intermediate broad Ha is 
about 700 light-days. Furthermore, the measured sizes 
based on observed total br oad Ha and H/3 are consi stent 
wi th the reported values in IKaspi et al.l ([2000, 2005J and 
in IPeterson et al.l (|2004ft . 

The clear size of intermediate broad component of Ha 
through CCF results indicate the intermediate BLR for 
Ha PG0052 can be mathematically determined. Based 
on the result, we try to discuss some unique characters 
based on intermediate BLR. 

3.5. Intermediate BLR 

Before to give some clear and further conclusion for 
mapped AGN PG0052. We discuss effects of intermedi- 
ate BLR on virial BH masses as follows. 

There are two parameters used for virial BH masses 
in Equation (1), line width and size of BLR. The two 
parameters are based on observational information. If 



intermediate broad component was accepted for Ha of 
PG0052, we will find that virial BH mass estimated 
from properties of total observed broad Ha (one inner 
broad component plus one intermediate broad compo- 
nent) should be some different from BH mass estimated 
from inner (or intermediate) broad component, and that 
estimated virial BH mass based on only inner (or inter- 
mediate) broad component of Ha should be more accu- 
rate than virial BH mass based on total observed broad 
Ha. Under the assumption of intermediate BLR ac- 
cepted for Ha of PG0052, we would find that 



J2. 

nHa 
-"■BLR 



fixR 



Hai 
BLR 



h x R 



Ha 9 
BLR 



(4) 



where <th 0i (^blr) an< 4 a Ha 2 (^blr) mean measured 
line widths (sizes of BLR) based on inner broad compo- 
nent and intermediate broad component of Ha, /i and 
ji represent parameter of flux weight: flux ratio of inner 
(intermediate) broad component to total observed broad 
Ha (/1 + /2 = 1). One simple but clear method to prove 
the equation above can be found in Appendix. 

It is very interesting that for PG0052 (only inner broad 
component for H/3 but inner broad plus intermediate 
broad component for Ha), virial BH masses based on 
total observed broad Ha should depends on properties 
of inner and intermediate component of Ha 



Mbh oc a\ a x R^ R 



x 



BLR/ 
2 

a Ha 2 



(5) 



x -^blrJ 



In the equation above, relation (j\ ai X i?BLR — c^ a2 x 

i?BLR. is accepted (BH mass estimated through proper- 
ties of inner broad Ha should be similar as mass through 
properties of intermediate broad Ha). Then we can com- 
pare BH masses based on observed total broad H/3 (inner 
broad component for H/3) and virial BH masses based on 
observed total broad Ha 



„2 v nHa 



BLR 



,R 



Hoi 
BLR 



R 



Ha2 

BLR 



a H(3 a -n B LR 



Vi + ffi + hM^ + ^gt) (6) 
it 



BLR 



R 



BLR 



In the equation above, the relation crj^g x i?BLR — 

CT Hai X-Rblr. — ff L 2 X-RgLii ^ s accepted (BH masses esti- 
mated through inner broad H/3 (i.e., observed broad H/3) 
should be similar as mass estimated on inner broad Ha). 
Thus, intermediate broad component of Ha for PG0052 
should lead to some different virial BH mass through to- 
tal observed broad H/3 (only one inner broad component 
in observed spectrum) and total observed broad Ha (one 
inner broad component plus one intermediate broad com- 
ponent), and lead to the result that PG0052 should be 

one outlier in the plane of o\aja\ a vs ^blr/^blr- 

Now let us check results of Equation (6) for PG 
0052+251. Flux ratio of inner broad Ha to intermediate 
broad Ha (total observed broad Ha minus intermediate 
broad Ha) is about 2.5 to 3 (2.5 from mean spectrum, 
and 3 for observed spectra), size ratio of intermediate 
BLR (for intermediate broad Ha) to inner BLR (for in- 
ner broad Ha) is about 6-7. Then the ratio of u\ a x -Rblr 



G 
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t° ^Hfl x ^blr should be about 2 based on Equation (6), 
which is well agree with position of PG 0052+251 in the 
Figure [T] In the Figure [TJ the area marked by dotted 
lines shows the accepted range for ratio of u\ a x R^r 
to cr^o x i?BL R under the size ratio of ^blr/^blr ^6 — 7, 
based on Equation (6). Furthermore, it is simple to check 
and confirm the correlations shown in Equation (4). The 
results indicate outlier PG 0052+251 in Figure [Jean be 
perfectly explained by properties of intermediate broad 
component of Ha, which further indicate the intermedi- 
ate broad component (or intermediate BLR) is reason- 
able. 

Before the end of the subsection, we can compare our 
measured size of intermediate BLR with reported results 
about size of intermediate BLR in literature. Zhu et al. 
(2009) have found one correlation between size of inter- 
mediate BLR and central BH masses (Figure 7 in Zhu ct 
al. 2009). If we accept M BH ~ (3.69 ± 0.76) x 10 8 M Q 
for PG 0052+251 (Peterson et al. 2004), then estimated 
size of intermediate BLR based on correlation shown in 
Zhu et al. (2009) is consistent with our measured value 
of ~ 10 18 cm, which provide further reliable evidence for 
intermediate BLR. 

3.6. Extended Part of Inner BLR? 

One question for existence of intermediate BLR is 
whether the component is only probable extended part of 
inner BLR, but not one isolated region. We answer the 
question as follows. Based on the results above, there 
are large distance between inner BLR and intermediate 
BLR, about 6001ight-days. If the probable intermediate 
BLR is just the extended part of inner BLR, the BLR 
of PG0052 should have a much extended size. So large 
extended size of BLR should smooth variations of ob- 
served broad emission lines, in other words, there should 
be no apparent variations of broad emission lines. We 
check the effects of extended size of BLR on variations 
by following mathematical procedure. 

Befor e proceeding fur ther, we accept the assumptions 
listed in [Peterson] ([1993D for reverberation mapping tech- 
nique, 1): continuum emission is from one central source 
which is much smaller than BLR, 2): both continuum 
emission and line emission are freely and isotopically 
propagating in central volume, 3): line emissions are in 
rapid response to ionizing continuum. Furthermore, we 
accept that line intensity from one region simply depends 
on number density of the region, / oc N(r) (some de- 
ta iled discussion for emission lines of AGN c an be found 
in lNetzer|[T990HOsterbroack fc Ferlandll2006[ ). Moreover, 
we simply accept that BLR is spherical shell (center of 
the sphere at position of central black hole) with depth 
about 600 light-days (inner layer radius is about 90 light- 
days, outer layer is about 700 light-days) . Then based on 
the geometrical structure, responsed light curve of emis- 
sion line based on input light curve of continuum e mission 
can b e created, as what we have done for 3C390.3 (| Zhang) 
120111) . The spherical shell (BLR) is firstly divided into M 
layers. As long as M is large enough, the effects of depth 
of each layer can be totally ignored. The line intensity 
from each layer (with radius r) can be determined by 



BLR, p ~ is accepted. If p is much large, the outer part 
of BLR should have few contributions to line emission. 
Thus, we select p = 0, there are similar flux strength from 
clouds in inner part of BLR and from clouds in outer part 
of BLR, which will show more clear effects of extended 
size of BLR on observed light curve of broad line. Once 
one layer meets the ionizing photos from central source, 
the line intensity of the layer is changed immediately as, 



Ii(t) oc h{t-l) x 



COUi (t) 
COUi(t — 1) 



(8) 



J(r) oc N(r) oc l/r p 



(7) 



where t is the date, Ii and coni mean line intensity from 
ith layer and arriving continuum emission (discrete data 
series with time separation of 1 day) for the layer. 

There is one point we should note that in our math- 
ematical procedure, the most simplest spherical shell 
geometry is assumed for BLR. Actually the geome- 
try is perhaps some different from tr ue struc t ure o f 
BLR of PG 0052+251 . As discussed inlGaskehl (T 2009); 
lEracleous efall ([2009f l: iSluse et al.l (|2011[ k iBentz et al.l 
(|2010a|) and reference therein, BLRs of AGN have a flat- 
tened distribution and that we always view them near 
pole-on, and BLR structures are very similar in most 
AGN. Although, the applied spherical shell geometry is 
oversimplified for PG 0052+251, results based on the 
simplified structure can still indicate true effects of ex- 
tended size of BLR on observed light-curve of broad line. 
Different structures of BLR should lead to some differ- 
ent overall trend of observed lig ht-curves of broad lines 
(such as one example shown in Zhang 2011), the large 
extended size of BLR still smooth observed light-curve 
of broad line. As one example to demonstrate effects 
of extended size of BLR, the oversimplified structure of 
BLR is significantly valid. 

Due to the simple procedure above, responsed light 
curve of Ha based on light curve of continuum emission 
of PG0052 can be created and shown in Figure [6] It is 
clear that the responsed light curve of Ha is very smooth, 
if extended size of BLR is about 600 light-days, i.e., if 
calculated intermediate BLR of PG0052 is just extended 
part of inner BLR. Furthermore, we also show one tested 
light curve of Ha, if BLR has smaller extended size about 
30 light-days. For the case with small extended size of 
BLR, some subtle features (arrows in the figure) shown 
in light curve of continuum emission can be reflected in 
light curve of emission line, however the features can not 
be found in created light curve of emission line for BLR 
with larger extended size. In the subsection, to clearly 
describe structures of BLR of PG0052 is not our objec- 
tive. Thus, we do not shown further results any more. 
The results shown in Figure M clearly indicate that the 
intermediate BLR can not be treated as extended part 
of inner BLR for PG0052, otherwise, the observed light 
curve of Ha should be much smooth. 

4. DISCUSSIONS 

After study of BLR of AGN for more than half one 
century, some information about geometrical structures 
of BLR has been mathematically determined by transfer 
function in the reverberation mapping technique through 
some special mat hematical methods (such as Maximum 
Entropy Method. iNaravan fc Nit vananda 1986; P eterson 



In order to find more clear effects of extended size of 1993; IPeterson et al . 1994; Home et al.lll991uGoad et al 
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1991 IWanders fe HorndH99l IPijpers fe Wanders! 1 1994 : 



KrolikI 119941: IWinge et all 119951: IBentz et all 1201011a : 
Sluse et al.ll2011l ). Certainly, besides the results through 
mathematical methods, there are some statistical re- 
sults about struct ures of BLR o f AGN (some reviews 
can be found in iGaskelll l2009t lEracleous et all 120091 : 
iDown et al.ll2010f) , based on properties of observed broad 
broad emission lines, such as the proposed model that 
one much broader component for wing of broad emis- 
sion line and one interm ediate broad component for core 
of broad emission l i ne (|Zhu et al.l 120091: IHu et al.l 2008; 
Sulentic et all 120001: iBrotherton et all 119941 iBrothertonl 
19961: iMason et al.lll996HBon et al.ll2009HPopovidl2007D . 
However, evidence is not so sufficient to confirm inter- 
mediate BLR of AGN. In this paper, besides the fit- 
ted results for broad Balmer emission lines, properties 
of line cores from PCA technique and the measured size 
of expected intermediate BLR are further confirmed ev- 
idence for intermediate BLR of PG 0052+251: besides 
one common inner component of BLR (inner BLR) with 
size about 100 light days, there is one other seriously 
obscured intermediate component of BLR (intermediate 
BLR) with size about 700 light days. 

Before proceeding further, we first simply consider one 
question what geometry is envisioned for intermediate 
BLR, where it suffers much more reddening than com- 
ponent giving rise to the inner BLR for PG 0052+251? 
Actually, it is not difficult to answer the question. Ob- 
scuration for intermediate BLR is not due to dust torus, 
but due to high density dust clouds (radial moving or 
not) between line-of-sight and intermediate BLR. The 
dust clouds have apparent effects on intermediate BLR, 
but no effects on inner BLR, due to small size of the 
clouds. The dust clouds can be confirmed by the spe- 
cial kind of AGN, AGN with their types changing be- 
tween type 1 and typ e 2, such as Mrkl018 dCohen et al. 
l9i aiGoodrid]H^90ll . NGC 76 03 (ITohline fe Osterbrock 



19791 . NGC 2622 and Mrk 609 (1Goodrichlll990D etc.. The 



large distance between inner BLR and intermediate BLR 
for PG0052 (about 6001ight-days) ensures enough space 
for the dust clouds. 

Furthermore , we consider the vitalization assumption 
(IGaskelll IlliM IWandel et all H99l IPeterson fe Wandell 
1999) for inte r mediat e BLR. As discussed in 
Brotherton et al.l (|1994l ) through simple photoion- 



ization model, properties of inner BLR and intermediate 
BLR are consistent with virialization assumption 



V 1 X i?BLR,l ~ V 2 X i?BLR,2 



(9) 



i?BLR means size of BLR (distance between BLR and 
central black hole) and V represents the rotating veloc- 
ity of broad emission line clouds in corresponding BLR 
(commonly, line width of broad emission lines, the sec- 
ond moment). For PG 0052+251, the sizes of two com- 
ponents of BLR and line widths of inner broad and inter- 
mediate broad components have been measured above, 



0"Ha 2 



4.7 + 0.4 



R 



Ho 2 
BLR 



R 



Hai 
BLR 



6 + 1.5 



(10) 



Here, the mean value (crHai/cHa 2 ) * s calculated by the 
measured line parameters for all observed spectra of PG 
0052+251, which is consistent with the one from mean 
spectrum of PG 0052+251. The value of RbZkI R mk. is 



calculated by the measured sizes of BLR listed in Table 
2. It is clear that the result is consistent with what we 
expect under the virialization method for BLR. 

Certainly, some effects on measured line parameters 
of mapped AGN s hould be discussed. As discussed in 
IKaspi et al.1 lf200 0) . the variable Fe II lines could al- 
ter line parameter measurements. However, Among the 
listed PG quasars in Table 1, PG 0052 and PG 0026 
are the only two objects without apparent optical Fell 
lines, the other 5 PG quasars have strong and appar- 
ent Fell emission lines. However, only PG 0052 is one 
outlier in Figure 1. Thus, effects of Fell on measured 
line parameters can be totally ignored. Another ques- 
tion we should discuss is siz e of BLR could ch ange, 
such as simple discussion in IKaspi et al.l (|2000l ) and 
detailed study a nd discussions for well-known mapped 
AGN NGC5548 (IBentz et al.l I2006L [20071 [Penney et al.1 
l2"Ol"0t IPeterson et al.l 120041: IWanders fe Petersonl 11999) 7 
However, we can find decreasing line width with increas- 
ing size of BLR, consistent with expected results under 
virialization assumption. Thus, effects of size of BLR 
changing could not explain the outliers in Figure 1. 

Before the end of the section, we should note that 
it should be interesting to find candidates for AGN 
with intermediate BLR by the outliers in the plane of 

-^blr/^blr vs ^h^/^Hq' especially for objects with 
much larger values of (crf^ x R^" R ) /{ay_a x -Rblr) ■ Thus, 
we will expect that mapped object NGC 4253 should be 
strong candidate with intermediate BLR. 



Finally, a simple summary is as follows. We first check 
properties of mapped AGN collected from literature in 
plane of cr^g/a^ vs ^blr/^-blr- Commonly, virial BH 
masses based on properties of observed broad Ha and 
H/3 should be coincident. However, among the mapped 
objects with measured sizes of BLR and line widths (sec- 
ond moment) based on long-period observed broad Ha 
and H/3, PG0052 and NGC4253 are two apparent out- 
liers in the plane of o"Hfl/ CT Ha vs ^blr/^blf.' which in- 
dicate BLRs of PG0052 and NGC4253 have some special 
characters. Then based on 55 public spectra of PG0052 
(jKaspi et al.ll2000h . the BLR of PG0052 is been carefully 
studied in detail. We find that line width ratio of total 
observed broad Ha to total observed broad H/3 is —0.7, 
which is much smaller than theoretica l /obser vational 
value of —0.9 found by iKorista fe Goadl (|2004f ) and in 
iGreene fe Hoi (|2005f) . Furthermore, the flux ratio of total 
broad Ha to total broad H/3 is about 6.8 (Balmer Decre- 
ment), which is not one reasonab le value (mean value of 
3.56 in iVanden Berk et al.ll200ll ) for BLUE quasar PG 
0052+251. Moreover, properties of line cores based on 
PCA technique indicate there is one inner broad com- 
ponent and one seriously obscured intermediate broad 
component in BLR. If the seriously obscured interme- 
diate BLR was accepted, properties of PG0052 in the 
plane of o\al<j\ a vs ^blr/^blr could be well repro- 
duced, which indicates that the intermediate BLR actu- 
ally is well appropriate to mapped quasar PG 0052+251. 
Finally, the large distance between inner component of 
BLR and intermediate component of BLR based on CCF 
results (about 6001ight-days) indicates the intermediate 
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component is not extended part of inner component of Dr. Kaspi S. to provide us the available observed spectra 
BLR. " of PG 0052 ( |http: / /wise-obs.tau ac.il/~shai/PG/| , and 
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APPENDIX 
TO PROVE EQUATION (4) 
Based on definition of second moment ((Peterson et al.ll2004f ) 

2 _ /A 2 xP A rfA f JXP x d\ \ 2 _ JX 2 xP x d\ 
° JP x dX [jl\dx) ~ JPxdX A ° (A1) 

where function P represents line profile, and A represents wavelength. For broad Ha of PG0052 which includes two 
components of Hai (inner broad component) and Ha 2 (intermediate broad component), we can find the correlation 
between line width of total broad Ha (cthq) and line widths of the two components of Ha (cjj ai and an a2 ) 

°Ha =/i x 4, + h x 4a 2 + A * ^(Ha x ) 

+ h x A 2 (Ha 2 ) - [fx x A (HaO + f 2 x A (Ha 2 )] 2 (A2) 
=/i x o\ ai + f 2 x cr| a2 + G[A (Hai), A (Ha 2 )] 

where f\ and f% are flux ratios of separated broad components of Ha to total broad Ha, i.e., /i = J Px,\dX/ J P\dX, 
/ 2 = / P\^dX/ J P x dX and /i + / 2 = 1. Ao(Hai) and Ao(Ha 2 ) are the first moments (center wavelengths) of the two 
components of broad Ha. It is clear that if inner component and intermediate broad component have not so much 
large different center wavelengths (Ao(Hai) ~ Ao(Ha 2 )), the equation a\ a ~ f\ X <7jj a + f-2 X cr^ a2 can be safely 
accepted. 

To mathematically prove equation ~ /i x + / 2 x i?^"^ should be very difficult. Here, we prove 

the equation by following Monte-Carlo method, based on completely HOMOGENEOUS light-curve s of continuum 
emiss ion (C(i)) and broad line emission (L(t)) of well-known reverberation-mapped AGN NGC5548 ([Peterson et alJ 
2002) collected from AGNWATCH project (http:// www.astronomy.ohio-state.edu/~agnwatch/). Light curves of inner 
and intermediate broad Ha are created through bootstrap method (one common used monte-carlo method to create 
a mock population from a given sample of data) applied to light curve of broad line H/3 of NGC5548. Observed light 
curve of NGC5548 including N data points can be described as [U, fluxi], where index i — 1,2, N means the ith 
data points in light curve, t and flux represent the observational date and corresponding flux density of emission line. 
In order to simply our mathematic procedure, the date U is re-created as isolated integral values with step of 1 day 
based on observed data series of NGC5548. Then mock light curves of inner and intermediate broad components of 
Ha can be created by the following two steps. On the first step, through the bootstrap method, a new sample of 
index from 1 to N, k = 1, ...N, is created. Certainly, there are some same values in sample of k. On the second step, 
mock light curve of inner broad Ha is created by [tj, fluxk], mock light curve of intermediate broad Ha is created 
by [ti + A, fluxk x f S ca], where A £ [2,160] days are integral values and represent distances between inner BLR 
and intermediate BLR, and f sca £ [0.2, 5] represents flux density ratio of intermediate broad Ha to inner broad Ha. 
Values of f sca ensure both the two broad components of Ha are apparent. Parameter A ensures that there is longer 
distance between intermediate BLR and central black hole than the distance between inner BLR and central black 
hole. Then the two steps above are repeated, until there are enough mock light curves. Furthermore, based on the 
bootstrap method, effects of simply different geometrical structures of inner and intermediate components of BLR 
have been s imply included (simple results about effects of geometric structures of BLR on measured size of BLR can 
be found in lZhu et al.ll200l^ . Based on the created mock light curves (date and flux density: [tinner, fluxk,inner] and 
[U, inter, f lu% k, inter]) of inner and intermediate component, it is easy to create the mock light curves ([ti^otjfluxk^tot]) 
of total broad line (inner component plus intermediate component) by 

^i,tot — ti.inter if ti.inter — ti, inner 

ti,tot Tejected if Winter 7^ ^i.inner (A3) 
fluX k ,tot = fluXk 

.inner + flux k ,tot if U 

.inter — ^i. inner 

fluXk iQt Vej GCted if ii,inter 7^ ^ynner 

Then values of -RblR' ^blr anc ^ -^blr can ^ e determined by CCF function. Here we notice that only CCF result 
with one apparent peak is accepted. Figure [7] shows the correlation (coefficient 0.97 with significance of deviation 
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TABLE 1 

Parameters for the 16 mapped AGN 



name 


"So 


^BLR 


<~H/3 


"BLR 


Public 




km/s 


light- days 


km/s 


light-days 




PG0026+129 


i f\s~* 1 I 1 n r~ 


98 - 1 -25.5 


1 7r7Q I O O F 

1773±28o 


11L0 -28.3 


Kaspi 


FGUUoz + ZOl 


1 OI Q _I_Q c 


163 - 7 - 38.3 


1 TOO 1 oc 

1 (00±OD 


on o + 24.5 

»y.»_ 24 A 


Kaspi 


t~> /~> no n i i 7£; 1 


on /i _L i qo 


183.6_ 13 3 


1 fl71 _l_ 1DK 

iy / i±iuo 


1 a a n+18.8 

146.9_ 18 9 


Kaspi 


IN Lj041o1 


O/l 00X70 


6.2_ 17 


1 OI A -LAO 
iyi4±4Z 


o. Iztl.o 


ALrlN WAl Url 


1~>/~^ 1/111 1 A A O 


2437±iyb 


n-i 7+36.0 
94 - 7 -31.5 


i /Tit I 1 an 

160 7±lb9 


1 o/i q+61.0 
124.3_ 61 7 


Kaspi 


FCt14zo+U15 


4z54zbz9U 


7 c c+30.5 
'3-0-32.5 


o /i I ono 

J442±d(J8 


nK n+29.9 
95.U_3 7 x 


Kaspi 


DPI 7 I 1 


o.ioo j_ i en 

z4oo±loU 


QA Q+19.1 

-25.2 


ZOZOltzil 


71 c-+29:6 
i _L . u 


Kaspi 


PP91 '30-1-000, 


1 491 4-P,0 


iys.4_ 23 4 




158 1+29- 8 


Kaspi 


NGC7469 


1164±68 


4 -'-1.3 


1456±207 


4 5+0.7 


AGNWATCH 


Mrkl42 


934±61 


2-90to.92 


859±102 


9 o 7 +0.76 
"*-°'_0.87 


NO 


SBS1116+583A 


1218± 9 f 


4 12+ 1 ' 41 


1528+.184 


9 oc+0.64 


NO 


Arpl51 


937±34 


a m +1.05 


1252±46 


4.081^0 


NO 


Mrkl310 


717+.75 


4 60+ ' 67 
4 - DU -0.62 


755+.138 


o 74 +0.60 


NO 


NGC4253 


726±35 


ZO.OU_Q.8g 


516±218 


R 94+1-65 
D -^-1.24 


NO 


NGC4748 


1035±74 


7 61 + 301 
'•Dl-4.64 


657±91 


5.63_ 2 25 


NO 


NGC6814 


1082±52 


M - oi -1.56 


1610+.108 


6 R7+0.88 
°- ' -0.90 


NO 



Notice:- 'Kaspi' means public observed spectra can be found from website 
http:/ /wise-obs. tau.ac.il/~shai/PG/ 

'AGN WATUH' mens public observed spectra can be found from the AGNWATCH project 
http: / /www. astronomy.ohio-state.edu/~agnwatch/ 

fMU' means there are no public observed spectra, the listed parameters are collected from Bentz et 
al. (2010). 

Size of BLR based on broad Balmer emission lines are the ones determined by the T cen t, and line 
widths are the ones measured from rms spectra as discussed in Peterson et al. (2004). 
We should note: listed values of parameters are values from more recent literature. 



TABLE 2 

Parameters of PG 0052+251 



Line 


a 


flux 


-Rblr(p) 


Rblr(c) 


RcCF(max) 




A 


10~ 16 crg/s/cm 2 


light- days 


light-days 




H/3 (broad) 


38.4 


1426.3 


109+.25 


111+.15 


0.76 


Ha (tot broad) 


37.4 


9612.5 


191±24 


214±19 


0.72 


Ha (inner broad) 


51.8 


5067.9 


110±20 


116+.16 


0.68 


Ha(int broad) 


24.1 


1649.5* 


703+.45 


678+.40 


0.38 



Notice:- Ha(tot broad) means observed total broad Ha. Ha(inner broad) means inner broad com- 
ponent of Ha. Ha(int broad) represents intermediate broad component of Ha. 

Second column gives the second moment of broad component in unit of A shown in mean spectrum of 
PG 0052+251, third column shows measured mean flux density of broad component for 55 observed 
spectra of PG 0052+251, fourth column shows measured size of BLR. through peak value of CCF 
results, fifth column shows the measured size of BLR through central value of CCF results, sixth 
column shows the maximum coefficient through CCF. 

Corresponding error of measured size of BLR is determined through bootstrap method as shown in 
Figure[5] 

1649.5* means the mean value is not so reliable and probably lower than the internal value, due to 
much weak intermediate broad Ha in some observed spectra. 



from zero -Pnuii ~ for 664 simulated data points) between i?gg R and x fx + R™^ x / 2 {f x = 1/(1 + / sca ) 

and f% = / S ca/(1 + /sca))- It is clear that the expected relation about size of BLR in Equation (4) can be commonly 
accepted. Furthermore, the results shown in Figure [7] indicate even there are two apparent two components in light 
curve of broad line, the CCF could be single-peaked not double-peaked appearance. 
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Fig. 1. — On the correlation between line width (second moment) ratio of ir^/ff^ and size ratio of ^blr/^blr' Standard error bar 
for data points is shown in bottom right corner. Object PG 0052+251 is marked with solid circle. The area marked by dotted lines around 
PG 0052+251 represents the estimated area in the plane, if intermediate BLR was accepted for PG0052 (results from equation (6) ). Solid 

line and dashed line represent the correlation ®\alv\ a = ^blr/^blr ant ^ confidence bands. For NGC4253, its true position 

([0.51, 4.1]) is much far away from the correlation. In order to clearly show the other objects, we use a upper arrow to represent that true 
position of NGC4253 should be much higher than the shown position [0.51, 2]. 




Fig. 2. — The mean spectrum of PG 0052+251, and best fitted results for emission lines around Ho and H/3. In top panel, solid line 
shows our mean spectrum created by PCA method (flux(5100A) = 1), dotted line shows the mean spectrum shown in Kaspi ct al. (2000) 
(/iua:(5100A) = 2). In bottom-left panel, dotted line shows the mean spectrum around H/3, thick solid line represents the best fitted 
results, dashed line represents broad component of H/3 fitted by one broad gaussian function, dot-dashed line shows narrow components of 
H/3 and [OIII] A4959, 5007A. In bottom-right panel, dotted line shows the mean spectrum around He*, thick solid line represents the best 
fitted results. Thick dashed line shows broad component, if one broad gaussian function is applied to fit Ha. Thin dashed line shows the 
inner broad component, dot-dashed line shows the intermediate broad component, if two broad gaussian functions are applied to fit Ha. 
Symbols © in the figure show positions of features of atmospheric A band near 7620Ain observed frame. 
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Fig. 3. — Properties of line cores. Solid line represents the first PCA component based on observed spectra with zero mean. Dashed line 
represents the mean spectrum based on observed spectra shown in Figure [2] Middle panel shows the comparison of line cores of Ha and 
H/3. In the panel, solid line represents line core of H/3 with flux density scaled by 3, dotted line represents line core of Ha. 
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Fig. 4. — Best fitted results for observed Ho. In each panel, dotted line shows observed spectrum, thick solid line represents the best 
fitted results, dashed line represents inner broad component of Ha, dot-dashed line shows intermediate broad component of Ha. In the 
case of Julian:2449596 (top-right panel), there is no apparent intermediate broad component. Symbols © in the figure show positions of 
features of atmospheric A band near 7620Ain observed frame. 
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Fig. 5. — Top-left panel shows results about Cross-Correlation Function (CCF). Thin solid line represents the result for broad 
H/3 (CCF(H/3, con • )). Dotted line represents the result for total broad Ha (CCF(Hat t,con °)). Thick solid line repre- 

scnts the result for inner broad Ha (CCF(Haj nner broad, con r 1 QQ^)) - Dot-dashed line shows the result for intermediate broad Ho 
(CCF(Ha 

int broadi con ° y) . Top right panel shows distributions of time-lags between total broad Balmer emission and AGN con- 
tinuum emission through bootstrap method, in order to compare our results with the results shown in Kaspi et al. (2000) and in Peterson 
et al. (2004). Bottom-left panel shows distributions of time-lag between inner broad Ha emission (after the subtraction of intermediate 
broad Ha) and AGN continuum emission through bootstrap method. Bottom-right panel shows distributions of time-lag between inter- 
mediate broad Ha emission and AGN continuum emission through bootstrap method. In the two bottom panels, solid line represents the 
distributions for time-lag determined by peak values of CCF results, dotted line represents the distributions of time-lag determined by 
central values of CCF results. 
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Fig. 6. — Effects of extended size of BLR on observed light curve of broad emission line. Solid line plus solid circle represents the observed 
light curve of continuum emission of PG0052, solid line represents the created light curve of emission line for the case that BLR with 
extended size about 30 light-days, dashed line represents the created light curve of emission line for the case that BLR with extended size 
about 600 light-days. The arrows show the positions of apparent features in light curve of continuum emission and in corresponding created 
light curve of emission line. 




Fig. 7. — On the correlation between RgL R and X /l + ^blh x f 2 ^ or t° ta l 664 simulated data points. Open circles represent 

444 data points with /2//1 > 1, and solid circles are for 220 data points with /2//1 < 1. Solid line represents the relation Rg" R = 



